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We analyze the applicability of our analytical theory of combination harmonics in canonical low- 
dimensional multi-band Fermi liquids, which was recently criticized by Champel (Phys. Rev. B 
65, 153403 (2002)). It is shown that his claim that our analytical theory does not apply at low 
temperatures and in clean samples, is incorrect. We demonstrate that the analytical theory of 
combination harmonics is in excellent agreement with the exact numerical results even at zero 
temperature and for clean systems, which are the most challenging for an analytical description. 
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The combination frequencies in the magnetization of 
the two-dimensional multiband metals with constant net 
electron density have been predicted, given simple phys- 
ical explanation, and described numerically in Ref. [1]. 
They were further numerically studied in Refs. [2,3] and 
observed experimentally by Shepherd et at [4] and by 
other groups [5]. We have also proposed an analyti- 
cal theory [6] of the combination Fourier components 
in the framework of the semiclassical Lifshitz-Kosevich 
approach [7]. More recently the theory has been ex- 
tended by taking into account the Dingle, spin and angle 
(Yamaji) reduction factors, and the nonquantized "back- 
ground" density of states [8] . The purpose of the present 
paper is to analyze our analytical theory in view of crit- 
ical remarks in Ref. [9], who claimed that "the chemical 
potential oscillations appearing in the arguments of the 
Fourier components were not taken into account? by the 
present authors. We attempt to clarify the relevant issues 
pertaining to our analytical expressions, in order to in- 
dicate explicitly the approximations made in the deriva- 
tions and resolve problems with their interpretation. Im- 
portantly, we also demonstrate that the analytical results 
for amplitudes of combination harmonics are numerically 
accurate even in most unfavorable circumstances at zero 
temperature and for clean samples. 

The basic equations of the theory [6] are those for the 
oscillating part F of the free energy, which is the ther- 
modynamic potential of the canonical ensemble, 

F = n( f ,)-(2p)- 1 (dh/dp) 2 , (1) 

and for the chemical potential p = po + p , where 
p = p^dQ/dfi [Eq. (12) and Eq. (9) of Ref. [6], respec- 
tively]. Here p = P~ X (N + J2 a P a ^ a ) ^ s * ne zero-field 
chemical potential, p and p a are total and partial zero- 
field densities of states (DOS), and A Q is the a— band 
edge. The oscillating part of the grand canonical poten- 
tial ft (Eq. (6) of [6]) is given by the standard expression 
[10] 



n( M ) = 2££(-irA;cosf^J, (2) 
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where A r a are the amplitudes of the single-band Fourier 
harmonics. For the sake of simplicity, we take the spin- 
splitting g— factors to be zero. Substituting this expres- 
sion into F, one obtains the combination amplitudes C r J a , 
as [5] 
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where to a — eB/m a . Our definition of the frequencies 
f a = 2irm a (p — A Q )/e in Eq. (2) and in Ref. [6] contains 
the exact chemical potential p rather than its zero-field 
value po. Hence, Eqs. (12) and (13) of Ref. [6] are exact. 
They fully take into account the chemical potential oscil- 
lations in the arguments of the Fourier components. 

Certainly, we did not consider our explicit expression 
for F, Eqs. (12) and (13) [6] as a final Fourier series, 
as should be obvious from our using the exact p, not its 
smooth part p$, in the expression for the free energy. The 
free energy F has been expanded in powers of p <C po as 

F^n(p Q ) + (2py 1 (dn(p Q )/dpo) 2 , (4) 

and, differentiating it with respect to magnetic field, we 
have obtained the result, Eq. (16) [6], for the magnetiza- 
tion amplitudes. The ratio of the combination M^ a , and 
single band (conventional) M\ Fourier amplitudes in a 
two-band metal was found at T = to be (Eq. (16) in 
[6]) 

*C = m a f a ± f a , 
M a m a + m a > f a 

It was generalized in Eqs. (21), (22) of Ref. [8] by taking 
into account a background DOS, the Dingle, spin, and 
Yamaji reduction factors. Obviously, the same result for 
the magnitude of this ratio can be obtained by differ- 
entiating our exact free energy, Eqs. (12,13) in [6], with 
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respect to the magnetic field B but keeping frequency f a 
constant. Therefore, the coefficients C r J a , in the second 
term of the free energy indeed yield the correct combi- 
nation frequency amplitudes in the magnetization, as we 
pointed in our original paper [6]. The expansion of the 
free energy in powers of ft described above is straight- 
forward, and it was not mentioned explicitly in [6]. The 
present discussion should clarify the point that the com- 
bination frequencies are fully defined by the second-order 
expansion coefficients C r J a , . 

Alternatively one can first differentiate the free energy, 
Eq. (1) with respect to B as [1], 

dF _(dh\ dhdfi ld& d on _ fon\ 
dB~\dBJ + l^~dB ~ ~p~d^'dB 1\l ~ [dB ) ' 

(6) 

which gives the exact expression for the magnetization 
(Eq. (5) of our original paper [1]). Then one can expand 
the result in powers of jl, as it was done by Champel in 
the second part of his paper [9] . Because the derivatives 
with respect to \x and with respect to B commute, the ap- 
proximate amplitudes should be the same as in Eq. (5). 
Indeed, the exact derivation of the magnetization, the 
approximate Fourier amplitudes, and the "main" result 
Eq. (14) of Ref. [9] are identical to our original expres- 
sions, Eq. (5) of Ref. [1] and Eq. (16) [6], Eqs. (21) and 
(22) [8], respectively. 

Champel mentions that "...the mechanism responsible 
for the combination frequencies in the Fourier spectrum 
of magnetization oscillations cannot apriori depend on 
the way the magnetization is derived, that is to say, on 
the use of a specific thermodynamic potential. Our fol- 
lowing goal is then to point out how the combination fre- 
quencies arise by considering directly the expression for 
the magnetization oscillations in the relevant thermody- 
namic limit" . Indeed, our results do not depend on the 
use of a particular thermodynamic potential. The mech- 
anism of novel combination frequencies [1] is based on 
the effect of chemical potential oscillations, which is im- 
portant in quasi-2D closed metallic systems and is absent 
in open systems. It can be derived with the use of either 
thermodynamic potential, its particular selection being a 
matter of convenience. To imply otherwise would be to 
misrepresent our work. Obviously, a natural choice of a 
thermodynamic potential for an open system (with a con- 
stant chemical potential) is (/i), while for a closed sys- 
tem with the constant number of particles TV this would 
be F(N). One is free to use fi(/x) at constant N as far as 
one accounts for the functional relation fi = fJ-(N) in the 
derivation, although it only makes the derivation cum- 
bersome without changing any results. The author of 
[9] seems to realize this, in fact using our formulas, e.g. 
Eqs. (3) and (4) in [9], to rederive our prior results. like 
e.g. (14) in Ref. [9]. 



The difference between the corresponding free energies 
Fand fl of the measured system is tiny, since it is propor- 
tional to the fluctuation of the carrier density, Eq. (1). 
However, the effect on susceptibility is greatly amplified, 
since two differentiations with respect to the field bring 
about very large factor (F/B) 2 » 1, Ref. [8], similar to 
the case of the usual de Haas- van Alphcn effect. Ob- 
viously, there will be no chemical potential oscillations 
when it is fixed by a reservoir, so that no combination 
frequencies due to these oscillations can be observed in 
an open (grand canonical) system. The reservoir com- 
monly implies a large system with an infinite continuum 
of nonquantized electron states. In this sense the role of 
the reservoir may be played by e.g. a one-dimensional 
(and, therefore, not closed) electron orbit with very large 
"background" density of states in the same quasi-2D sam- 
ple [2], observed in some cases. On the contrary, the un- 
usual combination frequencies appear when the system 
is closed, so that the number of carriers stays the same, 
and the chemical potential must oscillate [1]. In actual 
dHvA experiments the sample is usually measured while 
placed on non-conducting substrate with no electrodes 
attached, so the system is indeed closed. Note that the 
definition of the chemical potential does not require that 
the system be open. For normal Fermi liquids with the 
ground state energy Eo(k) for system of k particles the 
chemical potential [i = E {k + 1) — E (k) is uniquely 
defined for large k. 

Further, Champel writes about the origin of combi- 
nation frequencies: "... in 2D metals, the magnetization 
oscillations become more or less sensitive to the pres- 
ence of the chemical-potential oscillations and exhibit sig- 
nificantly different behaviors depending on the presence 
or absence of a finite reservoir of electrons. [2, 9, 10] In 
2D multiband metals this higher sensitivity to chemical- 
potential oscillations is expressed by the presence of the 
combination frequencies in the Fourier spectrum of mag- 
netization oscillations at very low temperatures, as shown 
numerically by Nakano [11]. Here, our aim is to prove 
analytically the existence of these combination frequen- 
cies in the magnetization oscillations" . One would note 
that the first claim is well known for single-band 2D met- 
als since the pioneering work by Peierls in 1933 [11]. One 
would easily see that the existence of combination fre- 
quencies has been proven numerically in our work [1] 
prior to Nakano's work. The analytical proof of their 
existence has already been given in Refs. [6,8]. 

Now we would like to address the claim of Ref. [9] 
that a the quantitative description of dHvA oscillations 
in terms of a Fourier series may break down and has to 
be done numerically" . Champel [9] has failed to mention 
that the numerical description of combination harmon- 
ics was done in our original [1] and subsequent papers 
[2,8]. To illustrate the accuracy of the analytical the- 
ory [5,7], let us analyze a simple two-band model. The 
results below show explicitly that the higher powers of 
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1. Relative ratio of the first approximation Mjj rst , 



Eq. (9), and exact M 11 combination amplitudes for a range 
of Dingle factors Rd- Note that the difference is in the range 
of 1 percent or less for both dirty and clean systems. We 
have assumed /2//1 = 7/11. Inset: the oscillating part of the 
chemical potential z = 2irjl/uj as a function of inverse mag- 
netic field 1/B. The two curves correspond to clean (the Din- 
gle factor Rd = 0.99, large amplitude of chemical potential 
oscillations 2nfl/uj) and dirty (Rd = 0.1, small amplitude) 
samples, respectively. 



fi -expansion can be neglected even at zero temperature 
and in a clean sample. Consider two bands with equal 
masses mi = m 2 = m and equal Dingle reduction fac- 
tors Rd = exp(— 2irT/ui), but with the different dHvA 
frequencies fx and fy. This two-band model is the most 
unfavorable case for the first-order expansion because in 
the presence of more bands the oscillations of the chemi- 
cal potential are reduced compared to the two-band case, 
as was correctly noted by Champel [9], The derivative 
of f2, Eq. (2) with respect to \x and the derivative of F, 
Eq. (6) with respect to B yield for T = 

00 fr> \r 

Z =-J2 - — [sin(rz + rfx/B) + sin(rz + rf 2 /B)] , 



r=l 



and 



M : 
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(8) 



x [/1 sin(rz + rfi/B) + f 2 sin(rz + rfy/B)} 



where z = 2tt/1/lu, and M is the oscillating part of the 
magnetization (to = eB/m). Expanding M in powers of 
z and neglecting z in the right hand side of Eq. (7) yields 
the amplitudes of the (/1 ± f%) Fourier components as [6] 
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(9) 



One can estimate the accuracy of this expression by 
Fourier transforming Eq. (8) without the expansion but 
leaving only the leading (r = 1) harmonics in the chemi- 
cal potential, z w Rd [sin(/i/B) + sin(/2/B)]. Then the 
ratio of the "exact" combination amplitude to the ap- 
proximate first-order amplitude is given by 
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where J n (x) is the Bessel function. This ratio is very 
close to unity at any value of the Dingle factor. In fact, 
the first order expansion [6] has even better accuracy 
than the estimate given by Eq. (10). 

For the numerical calculations we have used the ana- 
lytical expressions of the sums in Eqs. (7) and (8), 
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where Argiy stands for the argument of the complex num- 
ber w. The oscillating part of the chemical potential and 
the relative difference between the analytical and nu- 
merical amplitudes are shown in Fig. 1. The first-order 
Fourier amplitudes of the magnetization are equal to the 
numerical amplitudes within a few percent at any value 
of the collision broadening T. According to our theory, 
higher combination harmonics, like f\ + 2/2, should be 
exponentially suppressed compared to the leading har- 
monics /1 ± fi studied here. It is not surprising, there- 
fore, that the higher combination harmonics have not 
been clearly resolved even in very clean samples [12]. Evi- 
dently, the present results for the combination /i±/ 2 har- 
monics are sufficient to refute Champel's general claim 
that our analytical theory looses accuracy at low tem- 
peratures in clean systems. Indeed, the leading analyti- 
cal combination amplitudes [6,8] are very accurate even 
at low temperatures and in clean samples, contrary to 
claims in Ref. [9]. 

In conclusion, we have confirmed the validity of our 
analytical derivation [6], which does take into account 
the chemical potential oscillations. We have also shown 
that the experimentally observed analytical combination 
amplitudes [6] of the magnetization are numerically accu- 
rate even at zero temperature and even in clean samples. 
Therefore, the present analysis justifies the use of the 
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analytical theory [6,8] in the relevant range of all param- 
eters. 
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